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Abstract. We report STAR measurements of mid-rapidity yields for the A , A , K% , E~ , E + , f}~ , 
/) + particles in Cu+Cu and Au+Au ^/sjvjv = 200 GeV collisions. We show that at a given number of 
participating nucleons, bulk strangeness production is higher in Cu+Cu collisions compared to Au+Au 
collisions at the same center of mass energy, counter to predictions from the Canonical formalism. We com- 
pare both the Cu+Cu and Au+Au yields to AMPT and EPOS predictions, and find they reproduce key 
qualitative aspects of the data. Finally, we investigate other scaling parameters and find bulk strangeness 
production for both the measured data and theoretical predictions, scales better with the number partici- 
pants that undergo more than one collision. 
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1 Introduction 

Measurements of strangeness production in heavy-ion col- 
lisions were originally conceived to be the smoking gun of 
QGP formation [1] . It was argued that due to a drop in the 
strange quark's dynamical mass and increased production 
cross section, strangeness in the QGP would equilibrate on 
small time scales relative to a hadronic gas. Assuming a 
thermally equilibrated QGP hadronizes into a maximum 
entropy state, this hypothesis can be tested for heavy- 
ion collisions by comparing final state strangeness yields 
per participant to thermal model predictions from the 
Canonical formalism [2]. These predictions have qualita- 
tively reproduced various aspects of the data from Au+Au 
yj snn = 200 GeV collisions at RHIC, however, as with 
SPS energies, a complete theoretical description has yet to 
be achieved 3 . Furthermore, the strangeness saturation 
factor, 75, which characterizes the deviation in strangeness 
yields from thermal expectations, has been shown to ap- 
proach unity in central Au+Au ^/snn = 200 GeV colli- 
sions [4 . In these proceedings, we review the previously 
measured Au+Au ^/snn — 200 GeV data [4] [5] in the con- 
text of new data from approximately 40 million Cu+Cu 
\J snn — 200 GeV collisions recorded at the Relativistic 
Heavy Ion Collider in 2005. Measurements at AGS showed 
K + and K~ yields to be higher in lighter systems com- 
pared to the respective values in heavy systems at a given 
number of participants |6J, while measurements at the SPS 
showed higher K/ir ratios for the light systems also at a 
given number participants [7 . Whether these trends con- 
tinue up to RHIC energies, and what new information can 



be learned from strangeness production as QGP signature 
at RHIC, will be central issues for this analysis. 



2 Yield Extraction 
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The STAR Time Projection Chamber (TPC) is used to 
extract yields of A , A , K s , E~, S + , Q~ , and J? + as 
a function of transverse momentum, pt, via their dom- 
inant weak decay channels, which are A — > p + + 7r~, 

A -> p~ +7T+, Kg -> 7T++7T"", S — > A + TT~ , 5 — ► A + TT+, 

Q — ► A + K~ , and Q — > A + K + respectively. The de- 
cay products enter the TPC and are reconstructed using 
STAR's tracking software. The raw particle yields are cal- 
culated from the respective invariant mass calculations for 
the V0 and cascade candidates, and a combination of topo- 
logical, energy loss, and kinematic restrictions are placed 
to ensure the combinatorial background is minimal and 
can be described with a 2 nd order polynomial. To cal- 
culate the reconstruction efficiency, Monte Carlo particles 
are generated and propagated through a GEANT detector 
simulation. The generated charge clusters are then embed- 
ded into raw data, and the normal reconstruction process 
is applied. The efficiency is then defined as the ratio of 
reconstructed particles to the number of generated parti- 
cles, subject to the previously mentioned restrictions on 
the raw data. The efficiencies are applied to the raw data 
to give the corrected spectra. The A and A yields have 
contributions from the weak decay of charged and neutral 
S (anti) particles subtracted. Systematic uncertainty on 
the spectra points was found to be due to 1) An improper 
treatment of the vertex resolution in simulation, 2) Slight 
mismatch of the raw and simulated TPC hit distributions, 



2 



Anthony R. Timmins 



for the 



STAR Collaboration: Strange Particle Production at RHIC 



.v 



>■ Cu+Cu A 



r Au+Au A ' ^ 

,t< A 



- EPOS ver 1.67 

- AMPT ver 1.11 



10 



10^ 



<N„ 



10 



2 



6- 



4- 



• Cu+Cu S 

± Au+Au S 



_l I 



_l I I I I I Hi 



10 



1tf 



<M 



2 

0) 

>- 



2- 



• Cu+Cu7\ 
. ^ Au+Au ft 












■ 




■ 


10 


1 ° 2 <N part > 



• Cu+Cu" 
* Au+Au S 



_l I I I I 



10 



10 2 



2 

<a 
> 



3- 



q. 25 

Q. 
O 



■ -» Cu+Cu K° 




Preliminary 


Au+Au K° 




































■ 


10 




1 ° 2 <N > 

<|,, part > 



■J5 20 - Au+Au £2 + IT 



15 



10- 



" Cu+Cu £2 + 15 



I I 



_l I 



10 



10 2 



Fig. 1. Preliminary mid-rapidity per participant yields (dN/dy per (N pa rt)) of singly and multi strange particles for Cu+Cu 
and Au+Au collisions with ^/snn = 200 GeV divided by the respective p+p values. The A and A yields have been feed down 
subtracted in all cases. The green bars show the normalization uncertainties, and the uncertainties for the heavy-ion points are 
the combined statistical and systematic. The dotted and dashed lines show EPOS and AMPT predictions respectively where 
the p+p reference corresponds to the experimentally measured yields. Modulo the Q particles, the Cu+Cu trends (blue) in all 
cases lie above Au+Au (red) for a given particle species and model. 



and 3) Run-day variations in the yield due to variations in 
the TPC gain. Finally, we chose either an uit exponential 
or Maxwell-Boltzman in order to extrapolate to the low px 
region beyond the TPC acceptance so that dN/dy can be 
calculated for a given particle. More detailed descriptions 
of the yield extraction can be found here [8] [9] . 



3 Participant Scaled Yields 

Figure [1] shows participant scaled strangeness yields rela- 
tive to p+p values for Cu+Cu and Au+Au collisions with 
y / sjv r /v = 200 GeV. This variable is also known as the en- 
hancement factor as it characterizes the deviation from 
p+p collisions for participant scaled heavy-ion yields. The 
top panels show the singly strange particles which carry 
the majority of the strangeness produced, while the bot- 
tom panels show multi strange particles. For the singly 



strange particles it is clear that, for a given system, rises 
in participant yields are observed for each particle which is 
predicted by the Canonical formalism. However, at a given 
number participants above approximately 60 for Cu+Cu 
collisions, higher yields are observed. This is inconsistent 
with the Canonical formalism as it predicts a unique value 
for the enhancement factor at a given number of par- 
ticipants assuming the volume of the system is propor- 
tional to number of participants. Such a prediction also 
relies on a constant baryon chemical potential and con- 
stant chemical freeze out temperature for both systems 
which has been shown to be the case for Cu+Cu and 
Au+Au ^/sjvat = 200 GeV collisions within experimental 
uncertainty [lOj . Similar patterns are observed in the multi 
strange sector although any actual differences in Cu+Cu 
and Au+Au yields at a given number of participants are 
not as clear due to the larger fractional uncertainties for 
those measurements. 
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We also make comparisons to predictions from the 
AMPT [U] and EPOS [12] models which have been shown 
to describe the bulk features of hadron production well 
for Au+Au y^JV = 200 GeV collisions. The AMPT 
model is based on HIJING, and thus describes particle 
production in heavy-ion collisions via intra-nucleon string 
excitation and breaking (soft), and mini-jet fragmenta- 
tion (hard) where the excited nucleons fragment indepen- 
dently. For the settings used, mini-jet partons and newly 
produced hadrons can re-scatter, and without these mech- 
anisms, AMPT results reduce to that of HIJING. The 
EPOS model describes particle production with core and 
corona contributions. The core occupies the high partic- 
ipant density of collision zone (see figure [2]) and aims to 
mimic various QGP behavior. Once formed, it expands 
then hadronizes at a critical energy density similar to 
values predicted by Lattice QCD. The hadronization is 
treated via a statistical framework where strangeness is 
over-saturated with 7s = 1.3. This was chosen to the fit 
the Au+Au v^ivw = 200 GeV data by the EPOS au- 
thors. Corona production occurs in the low density region 
and can be thought of as a superposition of p+p collisions 
where strangeness is under-saturated at RHIC energies [S] . 
From peripheral (60-80%) Au+Au ^/JJ^n = 200 GeV col- 
lisions up to the most central, it was shown previously that 
core production is the major source of strangeness produc- 
tion where the relative contribution to total strangeness 
production rises from ~ 65% to ~ 100% [12]. Finally, for 
a given model, the default parameter set is used for all 
systems, centralities and particle species. For the singly 
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Fig. 2. Preliminary Monte Carlo Glauber calculation of par- 
ticipant densities (fm -2 ) in the transverse plane for central 
Cu+Cu collisions (left) and mid-central Au+Au collisions 
(right) with y / sjvjv = 200 GeV. The impact parameters for 
Cu+Cu and Au+Au are fixed to 2.37 and 9.27 fm respectively, 
and were chosen so both systems yield {N par t)—99. The radii 
at the redlines contain 99.7% of the respective Wood-Saxon 
distributions. 



strange particles, both the AMPT and EPOS models re- 
produce three key qualitative aspects the data: rises in 
per participant yields for a given system, higher yields at 
a given number of participants for mid-central and central 
Cu+Cu collisions compared to Au+Au, and a merging in 
per participant yields for peripheral Cu+Cu and very pe- 
ripheral Au+Au collisions. Quantitatively, EPOS is nearly 
always closer to the data than AMPT. With regard to the 
multi strange particles, both AMPT and EPOS also pre- 



dict the previously mentioned trends which are consistent 
with the data. EPOS is again closer to the Cu+Cu yields 
for all particles and over predicts the Au+Au yields, while 
AMPT under predicts both Cu+Cu and Au+Au yields. 

In figure [5] we show probability distributions for par- 
ticipant densities in the transverse plane for Cu+Cu and 
Au+Au yj s mn = 200 GeV collisions where the mean num- 
ber of participants are roughly the same in each case. The 
densities are clearly distributed differently owing to the 
different spatial geometries for the Cu and Au nucleons. 
In order to gain further insight into qualitative trends for 
strangeness production, the next section investigates the 
effects of geometry on the relations between other Glauber 
calculated quantities and strangeness yields for both the 
measured data and theoretical predictions. 



4 Empirical Scalings 

Given the similar systematic behaviour of all particles in 
figure [TJ we calculate the following quantity in order to 
concisely analyze strangeness production under various 
scaling hypotheses: 



dN, 



^ = 1.48^+4^*8 



dy 



dy 



dy 



dN A 
dy 



dNw 



dN~ 
, 4 — 

dy dy 

(1) 

which aims to approximate total strangeness production 
at mid-rapdity. The factors of 1.48 account for the pro- 
duction of (anti) S particles which have a similar quark 
content to the A particles, the factors of 4 accounts for the 
production of i*T , K Q Ll charged and neutral E particles 
[T3] . The fi particles are neglected due to incomplete cen- 
trality coverage in Au+Au ^snn = 200 GeV collisions, 
and they make a negligible contribution to the total s + s 
yield (less than 1% for central Au+Au). The left panel of 
figure [3] shows total mid-rapidity strangeness production 
per participant relative to p+p, where the relation in equa- 
tion [T] is applied to both measured data and theoretical 
predictions, as a function of the number of participants. 
As expected, total strangeness follows similar trends as 
the singly and multi strange particles. It is worth noting 
that while AMPT for the current settings always under 
estimates strange baryon yields, it nearly always over pre- 
dicts s + s (thus kaon production) yields which perhaps 
motivates improvements in the hadronization scheme. In 
particular, the introduction of baryon junctions has been 
shown to increase strange (anti) baryon production for 
other HI JING-based models [14] . 

Given the failure of the total strangeness production 
to scale with the number of participants, we test to see if 
mid-rapidity strangeness production scales with the num- 
ber binary collisions, (Nbi n ), therefore we plot dN( s+s ^/dy 
per (Npart) verses with v — 2(Nbi n ) / (N par t) in the mid- 
dle panel of figure [3] Such a scaling maybe indicative of 
large contributions from hard processes which is shown for 
charm quark production at RHIC [15j . However, although 
the Cu+Cu data moved slightly closer to Au+Au data, a 
common trend is neither observed for the measured data 
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Fig. 3. Preliminary mid-rapidity dN/dy for s + s per participant relative to p+p as a function of: the number of participants, 
{N par t)(ieit panel), v — 2(Nu n ) / (N par t) (middle panel) where Nti n is the number of binary collisions, and the fraction of par- 
ticipants that undergo more than one collision, / = {N par t>i) / {N par t) (right panel) where N par t>i is the number of participants 
that undergo more than one collision . 



nor the theoretical predictions. The middle panel of fig- 
ure [3] also provides a test for Kharzeev-Nardi scaling of 
strangeness production at RHIC. This decomposes soft 
and hard hadron production via the following relation: 



dN 
dy 



= (1 — x)n 



N, 



part 



VP 



pp-N bin 



(2) 



where n pp is the p+p yield, 1 — x is the fraction from soft 
processes (which scale with N part ), x is the fraction from 
hard processes (which scale with Nf,i n ), and < x < 1. It 
can be rearranged to give: 



dN -/N 



dy Impart 



PP 



/2 



x(v- 1) + 1 



(3) 



where the left term is equivalent to the value on the y-axis 
of figure [3J In the Kharzeev-Nardi framework, x should 
depend on the centre of mass energy (y / s]v r /v) in line with 
QCD predictions that the relative contribution of hard 
processes to particle production increases with increasing 
center of mass energy |16j . However, there is a collision 
system dependance (x appears larger for Cu) which is not 
naively expected in such a framework. 

Finally, in the right panel of figure [3] we test to see if 
strangeness production scales with the number of partic- 
ipants that undergo more than one collision, (N par t>i), 
therefore we plot dN( s+S )/dy per (N part ) verses / = 
(N P art>i) / (Npart) ■ For this scaling regime the measured 
data move closest together compared to the other scal- 
ing variables, while both EPOS and APMT Cu+Cu lines 



appear to lie on the respective Au+Au lines. The next 
section will discuss the potential relevance of the number 
of participants that undergo more than one collision for 
understanding strangeness production. 



5 Discussion 

Figure 3] shows the fraction of participants that undergo 
more than one collision as function of the number of par- 
ticipants for Cu+Cu and Au+Au ^snn = 200 GcV col- 
lisions. The distributions clearly resemble the measured 
data and theoretical predictions in figure [T] which is ex- 
pected from the scaling relations in the right panel of fig- 
ure [31 The success of this variable in the AMPT frame- 
work can be attributed to the fact participants with more 
than collision have more transverse momentum transfered 
to the valence quarks, and when these participants de- 
cay via string breaking, they produce more strangeness 
compared to participants with just one collision [T7]. As 
the fraction of these participants increases with central- 
ity for a given system, strangeness production per partici- 
pant also increases. When this fraction is higher in Cu+Cu 
compare Au+Au at a given number of participants, again 
strangeness production is higher for the former. 

Regarding the EPOS framework, the number of par- 
ticipants that undergo more than one collision is likely 
to have a direct relation to the size of core which, as 
mentioned, is the major source of strangeness production. 
Such a relation is expected as these participants are likely 
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to sit in the high density region of the collision as observed 
in figure[2l As for the fraction of participants than undergo 
more than one collision, it has been shown previously that 
core strangeness production per participant increases with 
the number of participants in Au+Au [T2]. Therefore a 
simple explanation for higher strangeness production in 
Cu+Cu at a given number of participants is a larger core 
compared to Au+Au. This would also be inferred from the 
higher number of participants that undergo more than one 
collision in Cu+Cu. 

Finally, we address how these observations possibly re- 
late to the original hypothesis of strangeness production 
as a QGP signature at RHIC. In this regard, the EPOS de- 
scription is closest as the core aims to represent deconhned 
matter with a unique hadron production scheme compared 
p+p. If strangeness equilibrates or indeed over-saturates 
in a QGP based on the original arguments, it is not incon- 
ceivable that 75 >= 1 upon hadronisation as mimicked by 
the core. On the other hand, the AMPT model offers an 
alternative view where soft strangeness production occurs 
in independently decaying nucleons rather than a coupled 
deconhned medium relevant to the original hypothesis. 
As observed in figure [3l this mechanism with the final 
state hadron interactions appears sufficient in producing 
the total number of mid-rapidity strange quarks. However, 
as also noted, the hadronization of these strange quarks 
leads to insufficient hyperon yield predictions compared 
to the measured data where EPOS generally does bet- 
ter. For further distinguishing power, it would interesting 
to see whether either model predicted the observed <j> en- 
hancement factor above one which is again contrary to 
expectations from the Canonical formalism at RHIC en- 
ergies [TS]- Indeed, an alternative core-corona approach 
has been shown to describe the <j> enhancement factors in 
Au+Au t/snn = 200 GeV collisions quite well [19] . 

6 Summary 

In these proceedings we have shown the systematic trends 
for mid-rapidity strange particle yields as a function of 
centrality for Cu+Cu and Au+Au ^/sjviv = 200 GeV col- 
lisions. We have found that both AMPT and EPOS mod- 
els reproduce the qualitative aspects for strangeness yields 
per participant with EPOS doing better quantitively for 
the strange baryon yields. It was also shown that despite 
the differing hadronization schemes, the successes of both 
models can be understood in the relatively simpler terms 
of nuclear geometry i.e. how the number of participants 
and the number of participants that undergo more than 
one collision evolve with centrality and system, which sug- 
gests this has a strong relevance to strangeness production 
at RHIC energies. 
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